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Abstract 


The  NAS A-Lewis  computer  program  for  the  calculation  of  com^plex  chemical  equilibrium 
compositions  (CEA)  is  proposed  as  a  means  of  characterizing  the  compositions  of  plasmas 
generated  via  an  ablating-capillary  arc.  Results  obtained  with  CEA  are  compared  to  those 
obtained  with  a  Saha  equation-based  (SE)  model  with  a  limited  reaction  product  set.  Case 
studies  conducted  at  temperatures  from  5,000  to  20,000  K  and  pressures  from  10  to  200  bar 
indicate  that  CEA  warrants  consideration  for  use  in  modeling  such  plasmas. 
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1.  Introduction 


Ablating-capillary  arcs  are  a  means  of  converting  electrical  energy  into  a  form  suitable  for 
the  ignition  of  propulsion  charges  in  large-caliber  gun  systems  [1].  In  ignition  systems  that 
utilize  this  approach,  electrical  energy  stored  in  a  pulse-forming  network  (PFN)  is  initially 
discharged  via  a  thin  metallic  wire  strung  through  the  capillary.  The  initial  discharge  explodes 
the  wire,  “instantaneously”  producing  an  ionized  (plasma)  gas  into  which  the  PFN  continues  to 
dump  energy.  Ablation  of  material  from  the  capillary  wall,  which  is  generally  attributed  to  (and 
modeled  as)  a  radiation  driven  process  [2],  serves  to  sustain  the  arc  by  replenishing  gas  that  exits 
the  capillary.  The  overall  process  produces  a  “dense”  (lO^^-lO^®  m'^),  “low-temperature” 
(10,000-30,000  K)  plasma  whose  composition  is  primarily  derived  from  decomposed  capillary 
material  [2].  This  fluid  subsequently  flows  into  the  propellant  bed  and  ignites  it  (directly). 

Ablating-capillary  arcs  were  one  of  the  first  electrothermal-chemical  (ETC)  ignition  concepts 
to  be  investigated,  and  experimentally  observed  attributes — including  short,  reproducible 
ignition  delays  [1],  propellant  bum  rate  enhancement  [3],  and  “temperature  compensation” 
[4,  5] — have  been  bases  for  believing  that  ETC-based  concepts  are  capable  of  substantially 
enhancing  the  ballistic  performance  (projectile  muzzle  velocity  or  kinetic  energy)  of  currently 
fielded  weapon  systems  [1].  Unfortunately,  the  mechanisms  that  underlie  the  experimental 
observations  have  yet  to  be  established,  making  the  design  of  a  fieldable  ignition  system  an 
elusive  goal.  This  report  summarizes  pan  of  an  effort  to  develop  a  state-of-the-art  computational 
fluid  dynamics  (CFD)  model  of  the  plasma-propellant  interactions  attending  ablating-capillary 
arc-based  ignition  [6,  7].  It  is  considered  that  such  a  model  will  serve  to  elucidate  the 
mechanisms  by  which  this  ignition  concept  yields  its  promise.  Moreover,  it  is  considered  that 
mechanisms  underlying  ablating-capillary  arc-based  ignition  will  be  relevant  to  (and  can  be 
exploited  for)  the  development  of  other  ETC-based  ignition  concepts. 

The  dynamics  of  the  generation  and  evolution  of  the  fluid  produced  by  an  ablating-capillary 
arc  are  highly  complex,  and  a  detailed  characterization  of  the  process  requires  cross-disciplinary 
expertise  in  plasma  physics,  fluid  dynamics,  and  chemical  kinetics.  The  work  presented  in  this 
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report  concerns  the  characterization  of  the  chemical  composition  of  the  effluent  from  an 
ablating-capillary  arc — an  issue  typically  addressed  by  reference  to  plasma  physics  concepts  and 
theories.  Of  particular  importance  to  the  characterization  of  an  ablating-capillary  arc  is  the 
concept  of  local  thermodynamic  equilibrium  (LTE)  (e.g.,  see  Griem  [8]).  The  concept  of  LTE  is 
relevant  to  systems  in  which  collision-induced  transitions  and  reactions  are  more  frequent  than 
radiative  transitions.  If,  in  addition,  there  is  microreversibility  between  the  collision  processes, 
the  steady-state  solution  of  the  rate  equations  yields  population  distributions  that  are  the  same  as 
those  of  a  system  in  complete  thermodynamic  equilibrium  (where  the  radiation  temperature  and 
the  temperature  of  the  reacting  particles  are  equal).  Thus,  when  LTE  can  be  assumed, 
Boltzmann  distributions,  Saha  equations,  and  the  mass  action  law  remain  valid  for  describing  the 
relative  populations  of  reacting  particles  even  though  the  radiation  distribution  may  deviate  from 
a  Planck  function.  These  formulae,  together  with  reaction  product  rovibronic  energy  level  data, 
provide  a  “microscopic”  description  of  the  plasma  that  allows  the  population  distributions  of  all 
species — neutrals  as  well  as  ions — to  be  calculated  based  on  knowledge  of  total  “heavy  nucleus” 
densities  and  temperature  (alone).  (A  model  based  on  these  formulae  is  subsequently  referred  to 
as  a  Saha  equation  [SE]  model.)  This  “reduction”  of  the  problem  greatly  facilitates  the 
characterization  of  a  plasma’s  chemical  composition. 

The  plasmas  produced  by  an  ablating  capillary  arc  are  relatively  dense,  making  the 
assumption  of  LTE  appropriate  for  such  systems,  and  SE  models  are  a  standard  means  for 
characterizing  their  chemical  compositions.  For  complex  systems,  the  main  effort  involved  in 
developing  an  SE  model  for  a  given  application  is  the  building  of  an  appropriate  reaction  product 
rovibronic  energy  level  database.  This  is  not  too  difficult  if  only  monatomic  species  need  to  be 
considered.  However,  at  the  low  end  of  temperatures  (10,000  K)  relevant  to  arc  discharge 
processes,  it  is  expected  that  the  density  of  polyatomic  species  will  be  significant  [9].  Since  the 
development  of  the  reaction  product  database  needed  to  model  such  systems  would  be  an 
extremely  arduous  and  time-consuming  task,  we  decided  to  determine  if  the  NASA-Lewis 
computer  program  for  the  calculation  of  complex  chemical  equilibrium  compositions  (CEA) 
could  be  employed  for  this  purpose.  CEA,  whose  origins  date  to  1967,  is  a  well-established  tool 
within  the  combustion  community  (where  it  is  often  simply  referred  to  as  the  NASA-Lewis 
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code).  Details  of  the  code  may  be  found  in  references  by  McBride  and  Gordon  [10,  11].  Briefly, 
the  program  permits  (among  other  things)  the  calculation  of  chemical  equilibrium  compositions 
for  assigned  thermodynamic  states.  It  does  so  based  on  a  “macroscopic”  mathematical 
approach — Gibb’s  free-energy  minimization.  One  of  the  attractions  of  CEA  is  the  extensive 
library  of  reaction  product  thermodynamic  data  (THERMO.INP)  that  has  been  developed  for  use 
with  it.  Information  on  more  than  1,300  species  is  included  in  the  library  in  our  possession. 
Originally  developed  for  characterizing  systems  at  temperatures  up  to  6,000  K,  in  recent  years 
the  library  has  been  expanded  to  allow  for  the  calculation  of  property  data  of  selected  species  at 
temperatures  up  to  20,000  K.  Among  the  species  with  properties  derivable  to  20,000  K  are  all  of 
the  neutral  and  +1  ions  of  the  atoms  and  diatoms  that  can  be  formed  in  a  gas  containing  C,  H,  N, 
and  O  nuclei.  This  expansion  of  the  library  raised  the  possibility  of  using  CEA  to  characterize 
ablating-capillary  arc  operation  in  this  temperature  regime,  and  thus  avoid  significant  model 
development  effort. 

With  respect  to  calculating  the  chemical  compositions  of  plasmas,  the  shortcoming  of  CEA  is 
that  it  doesn’t  include  a  correction  to  the  Gibb’s  free  energy  (calculation)  that  accounts  for  the 
influence  of  interactions  between  charged  particles.  These  interactions,  which  are  generically 
referred  to  in  this  report  as  “Coulomb  interactions,”  will  act  to  lower  the  ionization  potential  of 
plasma  constituents  and  (otherwise)  influence  constituent  thermodynamic  properties.  (That  is, 
the  thermodynamic  properties  of  a  plasma  constituent  will  depend  on  the  presence  and  density  of 
other  heavy-nuclei  in  the  plasma.)  Therefore,  it  is  not  clear  how  useful  CEA  will  be  in 
characterizing  the  composition  of  plasmas  generated  by  an  ablating-capillary  arc. 

To  determine  if  CEA  can  be  used  to  characterize  the  composition  of  the  effluent  from  an 
ablating-capillary  arc,  we  compare  the  compositions  computed  by  CEA  with  those  of  an  SE 
model  that  accounts  for  Coulomb  interactions,  but  employs  a  reaction  product  set  with  only 
monatomic  species.  The  SE  model  employed  in  this  study  is  the  same  as  that  employed  in 
Powell  and  Zielinski’s  (PZ)  model  of  an  ablating-capillary  arc  [2],  their  model  having  been 
employed  as  the  “source  term”  for  a  previous  CFD  model  of  the  expansion  of  a  plasma  jet  into  a 
combustion  chamber  [7].  The  comparison  is  conducted  for  plasmas  formed  from  the 
decomposition  of  polyethylene  [(C2H4)n],  a  material  that  has  been  extensively  used  and  studied 
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in  this  application  [2].  Results  are  presented  for  plasma  temperatures  in  the  range  from  5,000  to 


^V,WV  i.V  OiiU  ill  Ui^  lOilg^  liv^lli  J.V  LU  ^\J\J  uai. 


nn^ir*  r>Arvi'rvkA  4c  r>AlAT70T^f  i  rrrt  i +-1 

X  lllO  IWgllllW  lO  IWlWVCUlt^  U.XW  l^lXXLXV^Xl 


phase  of  an  interior  ballistic  cycle.  In  addition,  from  a  modeling  standpoint,  it  is  a  transitional 
regime,  with  CEA  expected  to  yield  better  estimates  for  lower  temperatures  and  pressures  and 
the  SE  model  expected  to  provide  better  estimates  for  higher  temperatures  and  pressures.  Thus, 
beyond  addressing  the  issue  of  employing  CEA  to  characterize  the  chemical  composition  of  an 
ablating  capillary  arc,  the  study  is  instmctive  as  to  the  use  of  an  SE  model  for  characterizing  low- 
temperature  plasmas  when  the  reaction  product  data  set  employed  with  the  model  includes  only 
monatomic  species. 


2.  Model  Considerations 

2.1  The  SE  Model.  For  a  plasma  in  LTE,  the  relative  concentrations  of  the  (J+lf'  and  f' 
ionization  states  of  the  i*  atomic  nucleus  at  a  temperature  [T)  may  be  estimated  via  the  Saha 
roi 
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where  nij  is  the  number  density  of  specie  is  its  “isolated  state”  ionization  potential,  and 

A  E“j  is  the  reduction  of  E,"  due  to  Coulomb  interactions.  Other  parameters  in  the  equation  are 

the  electron  number  density  (rig),  the  partition  function  [ZiJT)],  the  rest  mass  of  an  electron  (mg), 
Planck’s  constant  (h),  and  Boltzmann’s  constant  (kb).  For  monatomic  species,  the  partition 
function  is  a  function  of  a  specie’s  electronic  energy  level  structure  only  and  is  computed  using 

Z.j(T)  =  ^g.j^,exip 
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where  gij^k  is  the  statistical  weight  (factor)  for  energy  level  EZ  . 
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In  applying  equation  (1)  for  the  computation  of  chemical  compositions,  the  estimations  of 
AE^j  and  ZiJT)  are  issues.  The  parameter  AE“j  reflects  the  fact  that  Coulomb  interactions 

allow  the  electrons  of  an  atom  or  ion  to  escape  the  central  force  field  of  the  nucleus  at  energies 
less  than  the  ionization  potential  of  an  isolated  atom  or  ion.  Thus,  the  ratio  of  (j+1)*^  to 
ionization  states  for  the  nucleus  in  a  plasma  will  be  higher  than  that  estimated  assuming  that 
charged  constituents  do  not  “interact.”  For  estimating  AE“j,  we  chose  to  follow  PZ  [2]  and 
calculate  this  parameter  based  on  the  model  proposed  by  Eberling  and  Sandig  [12], 


AE°°.  = - ^ - 
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where  e  is  the  electronic  charge,  is  the  permittivity  of  free  space,  A  is  the  DeBroglie 
wavelength. 
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and  Ad  is  the  Debye  length. 
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The  parameter  z  in  equation  (5)  represents  the  effective  charge  on  an  ion,  this  parameter  being 
computed  via  the  equation 


z  = 
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where  xij  is  the  ratio 
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It  is  worth  noting  that,  in  this  approach,  increases  nearly  in  proportion  to  while 

decreasing  in  proportion  to  T^'^. 

The  difficulty  of  computing  can  be  appreciated  by  noting  that  the  partition  function 

becomes  infinite  unless  the  summation  is  truncated.  (Near  the  series  limit,  the  exponential  term 
remains  nearly  constant  while  the  statistical  weight  increases  with  principle  quantum  number 
[«].)  To  address  this  difficulty,  it  is  typically  assumed  that  Coulomb  interactions  have  a 
negligible  effect  on  the  energy  level  structure  of  individual  constituents  and  the  reduction  in 
ionization  potential  is  employed  as  a  basis  for  limiting  the  summation  in  equation  (2)  to  energy 
levels  [8] 


E’^  kE'  -  AE-".  (8) 

With  this  set  of  assumptions,  calculation  of  Z,/!)  simply  becomes  a  matter  of  identifying  all  of 
the  energy  levels  for  a  constituent  up  to  this  limit.  Identification  of  these  levels  can  be 
accomplished  by  reference  to  experimentally  based  data  tables  and/or  through  computational 
techniques.  We  have  relied  on  both  approaches,  the  details  of  which  are  given  in  Appendix  A. 

Estimates  of  chemical  composition  follow  from  the  solution  of  equations  (l)-(8),  which  are 
identical  to  those  employed  by  the  PZ  model  for  simulating  an  ablating-polyethylene  capillary 
arc  [2].  For  the  calculation  of  the  partition  functions  in  this  study,  however,  we  include  higher 
energy  levels  for  C,  C"^,  and  H.  In  the  case  of  carbon — the  only  case  where  the  inclusion  of 
higher  energy  levels  makes  a  significant  difference — ^PZ’s  study  was  based  on  calculations  that 
included  levels  with  energies  up  to  71,500  cm  *.  This  proves  to  be  adequate  for  their  study 
because,  through  equation  (8),  the  electron  number  densities  generated  by  the  arc  effectively 
eliminate  the  need  to  consider  higher  energy  levels.  However,  partition  functions  computed  with 


this  cutoff  begin  to  diverge  significantly  from  values  tabulated  by  Drawin  and  Felenbok  [13]  for 
temperatures  above  about  15,000  K  unless  exceeds  about  4,000  cm  *  (0.5  eV).  [From 

equation  (4),  ««  must  exceed  about  4xl0^'‘m'^  at  15,000  K  to  produce  a  AE'°j  of  this 

magnitude.]  A  AE,”  of  4,000  cm'*  truncates  levels  with  energies  greater  than  86,9(K)  cm'*,  and 

the  partition  function  calculations  for  carbon  prove  to  be  relatively  insensitive  to  the  inclusion  or 
omission  of  its  levels  that  exist  between  71,500  and  86,900  cm'*.  [Energy  levels  with  principle 
quantum  numbers  equal  to  3,  4,  and  5  exist  in  this  range.] 

In  the  interest  of  avoiding  database-related  limitations  beyond  those  we  specifically  intended 
to  impose,  more  wide-ranging  agreement  with  Drawin  and  Felenbok  [13]  was  sought. 
(Unfortunately,  Drawin  and  Felenbok  do  not  provide  sufficient  detail  about  the  bases  for  their 
calculations  to  do  them  directly.)  To  get  better  agreement,  energy  level  data  for  carbon  that 
included  all  predicted  levels  up  to  n  =  10  were  employed.  This  information  was  primarily 
obtained  from  references  by  Moore  [14]  and  Striganov  and  Sventitskii  [15],  but  also  includes 
predicted  levels  not  found  in  either  of  these  works.  These  missing  levels  were  added  based  on  a 
“fill”  procedure  that  is  described  in  Appendix  A.  Partition  functions  computed  from  Table  A-1 
(in  Appendix  A)  are  in  good  agreement  with  those  tabulated  by  Drawin  and  Felenbok  for  all 
AE,"  and  temperatures  from  2,750  to  20,500  K. 

We  also  note  that,  as  expressed  in  equations  (l)-(8),  the  concentrations  of  e',  C,  C'"'',  H, 
and  !!*■  are  functions  of  the  temperature  and  electron  number  density  (r.n^),  not  temperature  and 
pressure  {T,P)  per  se.  To  determine  concentrations  as  a  function  of  temperature  and  pressure, 
and  thus  facilitate  direct  comparison  between  the  SE  model  and  CEA,  we  developed  an  iterative 
procedure  to  find  the  electron  number  density  at  a  desired  temperature  and  pressure.  The  first 
step  in  the  procedure  involved  setting  z  equal  to  1  (based  on  the  expectation  that  the 
concentration  of  C"*'  would  be  extremely  small)  and  calculating  AE;*  and  Zi/T)  for  two  values 

of  Tie  that  bounded  the  electron  density  at  the  desired  pressure.  For  each  (n^,  z)  pair,  species 
concentrations  were  derived  and  an  updated  estimate  for  z  obtained.  The  new  {rig,  z)  pair,  in  turn, 
led  to  a  new  estimate  for  the  chemical  composition  (for  the  given  electron  number  density). 


7 


Succeeding  estimates  were  obtained  until  the  chemical  composition  and  z  were  consistent  to  nine 
significant  digits. 

Once  the  chemical  compositions  of  the  bounding  (n^,  z)  pairs  were  established,  the  pressures 
corresponding  to  these  pairs  were  found  from  the  ideal  gas  law, 


P  =  kTY, 


(9) 


A  bisection  root-finding  technique  was  then  employed  to  obtain  a  new  estimate  for  ««  (and 
chemical  composition)  whose  pressure  was  closer  to  the  desired  value.  The  convergence 
criterion  for  the  root  finding  routine,  which  is  slightly  modified  version  of  a  program  given  by 
Press  et  al.  [16],  was  set  at  1  x  10^^  m'^.  This  criterion  was  selected  based  on  preliminary  work 
that  showed  that  the  electron  number  densities  of  plasmas  considered  in  this  study  would  be 
greater  than  7  x  10^^  m"^.  Moreover,  when  this  criterion  is  met,  the  pressure  of  the  composition 
is  found  to  agree  with  the  desired  pressure  to  better  than  six  significant  digits. 

2.2  CEA.  The  most  recent  version  of  CEA  of  which  we  are  aware  (i.e.,  CEA600)  was 
released  in  1998,  but  we  had  trouble  getting  it  to  run  a  “TP”  calculation  (i.e.,  calculate  a 
composition  for  an  assigned  temperature  and  pressure  state).  Therefore,  we  reverted  to  CEA300, 
which  was  released  in  1997.  Like  the  SE  model,  CEA  employs  the  ideal  gas  law  to  model  the 
equation  of  state.  For  the  purposes  of  this  study,  we  modified  CEA  so  that  it  could  be  called  as  a 
subroutine  by  a  program  that  handled  input  and  output  chores.  We  also  added  thermodynamic 
property  data  (coefficients)  for  to  THERMO.INP.  Details  of  the  inclusion  of  this  specie  to 
the  data  file  are  described  in  Appendix  B. 


23  Temperature  and  Pressure  Bounds.  The  two  models  are  compared  for  temperatures 
ranging  from  5,000  K  to  20,000  K  and  pressures  from  10  bar  to  200  bar.  The  low  end  of  the 
temperature  range  was  set  in  reeognition  of  the  fact  that  CEA  is  well  established  for  temperatures 
up  to  6,000  K.  At  such  temperatures,  the  concentration  of  charged  particles  in  a  gas  composed 
from  C,  H,  N,  and  O  nuclei  will  be  very  small  and  Coulomb  interactions  will  be  negligible. 
Thus,  it  is  expected  that  CEA  will  provide  good  estimates  at  such  temperatures.  The  high  end  of 
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the  temperature  range  was  based  on  the  limit  of  applicability  of  the  data  in  THERMO.INP.  The 
pressure  range  is  applicable  to  ablating-capillary  arc  dynamics  during  the  ignition  phase  of  the 
interior  ballistics  cycle.  Higher  pressures,  which  would  be  pertinent  to  ETC  propulsion  concepts 
in  which  the  plasma  is  injected  later  in  the  interior  ballistics  cycle,  were  not  considered  because 
the  ideal  gas  equation  of  state  underlying  both  models  becomes  a  questionable  assumption.  The 
calculation  of  chemical  concentrations  at  high  pressures  requires  the  use  of  a  program  (such  as 
BLAKE  [17])  that  is  designed  to  address  nonideal  behavior.  The  influence  of  the  equation  of 
state  on  the  chemical  composition  of  an  arc  is  beyond  the  scope  of  the  present  work. 

3.  Results  and  Discussion 

As  developed  to  date,  CEA  and  the  SE  model  have  different  inherent  limitations:  (1)  CEA 
does  not  have  an  ability  to  account  for  Coulomb  interactions,  and  (2)  the  SE  model  has  a  reaction 
product  data  set  that  is  limited  to  monatomic  species.  Moreover,  they  employ  different 
mathematical  approaches  for  establishing  equilibrium.  Coupled  with  the  fact  that  CEA  is 
virtually  a  “black  box,”  these  differences  have  the  potential  to  obscure  bases  for  differences  in 
their  results.  To  address  this  issue,  we  first  compare  results  from  the  models  under  a  set  of 
imposed  limitations  that  highlight  “real”  differences  between  them.  As  points  of  reference  for 
the  comparisons  to  be  shown,  model  computations  for  the  temperature  dependence  of  plasma 
compositions  at  a  pressure  of  50  bar  are  shown  in  Figures  1  and  2. 

As  a  first  comparison  of  the  two  approaches,  the  reaction  product  set  employed  in  CEA  is 
limited  to  those  species  (e',  C,  C"^,  C"^"^,  H,  and  H"^)  considered  in  the  SE  model,  and  the  ionization 
potential  reduction  parameter  {A  E“j )  of  the  SE  model  is  set  equal  to  0.  Modified  in  this  manner, 

the  two  approaches  (should,  in  principle,)  become  mathematically  equivalent.  As  such,  it  was 
expected  that  the  models  would  yield  the  same  results.  As  shown  in  Figure  3,  this  expectation  is 
(essentially)  realized  for  all  species  except  for  temperatures  from  5,000  to  15,000  K,  the 
difference  in  predicted  species  concentrations  being  less  than  5%  over  this  temperature  range. 
(In  view  of  the  fact  that  C"^"^  concentrations  are  more  than  3  orders  of  magnitude  less  than  any 
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Figure  1.  SE  Model  Calculations  of  Reaction  Product  Number  Densities  vs. 

Temperature  for  a  50-bar  Gas  Formed  From  Decomposed  Polyethylene. 
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Figure  2.  CEA  Calculations  of  Reaction  Product  Number  Densities  vs. 

Temperature  for  a  50-bar  Gas  Formed  From  Decomposed  Polyethylene. 
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other  specie,  the  results  for  can  be  discounted.)  However,  as  the  temperature  increases 
above  this  range,  the  predictions  of  the  two  models  begin  to  diverge  significantly. 

We  suspected  that  the  differences  observed  as  the  temperature  increases  above  15,000  K  are 
due  to  differences  in  the  number  of  energy  levels  employed  to  compute  the  thermodynamic 
property  data  of  carbon  given  in  THERMO.INP  and  those  employed  to  compute  carbon’s 
partition  function  for  the  SE  model.  This  suspicion  is  borne  out  by  the  results  shown  in  Figure  4, 
which  compares  (1)  carbon’s  heat  capacity  as  computed  by  the  THERMO.INP  coefficients,  and 
(2)  the  approach  outlined  in  Appendix  B,  where  energy  levels  with  prineiple  quantum  numbers 
up  to  n  =  10  have  been  used  for  the  computation.  To  verify  that  the  difference  is  due  to  fewer 
energy  levels  being  employed  in  the  procedure  to  generate  the  coefficients  in  THERMO.INP, 
heat  capacity  values  calculated  when  energy  levels  above  71,500  cm''  are  excluded  are  shown 
for  comparison.  This  result  suggests  that  the  “TEMPER  method”  [18]  has  been  employed  to 
derive  the  THERMO.INP  coefficients  for  carbon  atoms.  (In  this  method,  E“j  is  lowered  as  a 

function  of  the  temperature.)  If  the  limited  energy  level  set  employed  to  generate  Figure  4  (c)  is 
employed  for  SE  model  calculations  of  chemical  composition,  the  results  obtained  from  the  two 
models  are  in  good  agreement  for  the  entire  temperature  range  considered  (see  Figure  5). 
(Again,  the  relatively  large  differences  observed  for  C"’"’’  are  related  to  the  fact  that  it  is  present  in 
trace  amounts.) 

The  results  shown  in  Figure  5  suggest  the  equivalence  of  truncating  the  partition  function 
summation  in  the  SE  model  [per  equation  (8)]  and  restricting  the  number  of  levels  used  for  the 
calculation  of  thermodynamic  properties.  The  cutoff  expressed  in  the  thermodynamic  property 
data  does  not,  however,  lower  (or  C'^'’’)  energy  levels  relative  to  C.  That  is,  CEA  does  not 
have  a  counterpart  to  the  influence  of  A  E^j  in  the  SE  model  as  realized  through  its  inclusion  in 

the  exponential  term  of  equation  (1).  This  limits  the  correspondence  that  can  be  achieved 
between  the  two  approaches  by  simply  limiting  the  energy  levels  considered  in  the  generation  of 
the  thermodynamics  data  contained  in  THERMO.INP.  This  lack  of  correspondence  is 
observable  in  the  results  presented  in  Figure  6,  which  compares  results  obtained  when  the 
reaction  product  set  employed  by  CEA  is  limited  to  those  species  considered  in  the  SE  model 
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Figure  4.  The  Dimensionless  Heat  Capacity  of  Neutral  Carbon  Atoms  as 
Computed  by  (a)  the  Method  Outlined  in  Appendix  B  With 
Energy  Levels  up  to  89,800  cm'*  Included,  (b)  the  Coefficients  in 
THERMO.INP,  and  (c)  the  Method  Outlined  in  Appendix  B  With 
Energy  Levels  up  to  71,500  cm  *  Included. 
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and  the  SE  model  calculations  are  performed  as  done  for  Figure  5,  but  with  equation  (4)  now 
“turned  on.”  The  comparison  shows  that,  as  expected,  CEA  predicts  lower  concentrations  of 
and  than  does  the  SE  model.  This,  in  turn,  influences  the  results  obtained  for  e-,  H,  and  H"^. 

Figures  7-9  compare  results  from  the  two  models  when  no  limitations  are  imposed  (beyond 
those  inherent  in  their  development  to  date).  These  comparisons  show  that  the  models  predict 
significantly  different  concentrations  of  monatomic  species  for  temperatures  below  about 
8,000  K.  (The  point  at  which  they  begin  to  disagree  depends  on  the  pressure.)  This  difference  is 
attributable  to  significant  concentrations  of  polyatomic  species  being  predicted  by  CEA  in  this 
temperature  range  (see  Figure  10).  Given  that  both  models  predict  that  electron  number  densities 
for  temperatures  less  than  10,000  K  will  be  less  than  6x  10  m'  ,  AFT  will  be  less  than 

0.15  eV.  Thus,  Coulomb  interaetions  will  have  a  negligible  influence  over  the  equilibrium.  In 
view  of  these  considerations,  we  expect  CEA  to  provide  estimates  that  are  superior  to  those  from 
the  SE  model  in  the  temperature  range  below  8,000  K. 

4.  Considerations  Related  to  Modeling  Ablating-Capillary 
Arc  Discharge  Processes 

The  results  presented  highlight  the  conundrum  facing  a  modeler  who  wishes  to  use  CEA. 
The  biggest  asset  of  CEA  is  the  large  reaction  product  database  that  has  been  developed  for  use 
with  it.  However,  the  bases  for  the  data  for  temperatures  greater  than  6,000  K  are  relatively 
inaccessible,  and  it  is  more  work  to  check  the  data  than  it  is  to  generate  it.  Thus,  CEA  only 
reduces  model  development  effort  if  the  modeler  is  willing  to  simply  accept  the  data  in 
THERMO. INF.  Given  the  sensitivity  of  the  results  to  the  method  employed  to  generate  the 
THERMO. INP  data,  such  an  approach  makes  us  uncomfortable.  However,  from  the  standpoint 
of  developing  a  detailed  model  of  the  evolution  of  the  plasma  from  an  ablating  capillary  arc,  the 
results  presented  in  Figures  7-9  provide  a  measure  of  justification  for  using  CEA  to  calculate  the 
chemical  compositions  of  such  plasmas.  Moreover,  they  indicate  the  limitations  of  an  SE  model 
that  utilizes  a  reaction  product  database  that  includes  only  monatomic  species. 
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Figure  10.  CEA  Calculation  of  the  Fraction  of  Total  Heavy  Nuclei-Containing 
Particles  That  Are  Monatomic. 
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As  an  initial  step  toward  developing  a  CFD  model  of  ablating-capillary  arc  ignition,  we  are 
using  the  PZ  model  of  arc  dynamics  to  calculate  the  temperature  and  pressure  (state)  of  the 
effluent  at  the  capillary  exit  plane  and  using  CEA  to  calculate  the  composition  for  this  state.  The 
PZ  model  calculations  of  constituent  densities  are  disregarded.  Of  course,  the  dynamics  of  the 
ablating  capillary  arc  process  depend  on  the  constituent  densities  in  the  arc,  and  thus  this 
approach  is  internally  inconsistent.  An  aspect  of  this  inconsistency  is  shown  in  Figure  11,  which 
compares  the  effective  molecular  weight  of  the  plasma  at  the  same  temperature  and  pressure  as 
computed  by  CEA  and  the  SE  model.  However,  this  inconsistency  is  considered  tolerable  when 
judged  from  the  standpoint  of  developing  a  detailed  model  of  the  discharge  process.  By  its 
nature,  such  a  development  effort  needs  to  be  done  incrementally,  and  CEA  has  allowed  us  to 
make  timely  preliminary  assessments  of  the  impact  of  polyatomic  species  on  the  evolution  of  the 
plasma  outside  the  capillary  [6].  The  results  obtained  with  this  approach  are  also  providing 
guidance  for  efforts  to  build  the  reaction  mechanisms  and  transport  coefficient  data  needed  for 
the  CFD  model  [6,  19].  And  such  benefits  will  have  even  greater  importance  when  systems  with 
N-  and  O-containing  species  are  considered,  their  chemical  compositions  being  much  more 
complex  than  the  “hydrocarbon”  plasmas  considered  here. 

In  the  course  of  this  study,  an  approach  for  incorporating  the  effect  of  Coulomb  interactions 
into  a  Gibb’s  free-energy  minimization  calculation  of  chemical  compositions  was  identified 
[20,  21],  and  we  consider  that  CEA  could  be  modified  to  incorporate  this  approach.  However, 
the  data  in  THERMO.INP  would  probably  also  need  to  be  modified,  effectively  voiding  the 
potential  benefit  of  such  an  undertaking.  Moreover,  CEA  is  relatively  inconvenient  to  use  as  a 
subroutine,  and  trying  to  modify  the  PZ  code  to  utilize  it  (and  thus  resolving  the  inconsistency)  is 
not  an  attractive  option.  Therefore,  as  the  next  step  in  code  development,  we  intend  to  develop  a 
new  model  of  ablating  capillary  arc  dynamics  that  is  integral  to  the  CFD  model  and  to  write  our 
own  subroutine  for  chemical  composition  calculations.  (The  mathematical  approach  we  will 
employ  remains  to  be  decided.)  This  effort  will  resolve  the  inconsistency  inherent  in  the  current 
approach  [6],  and  it  will  facilitate  our  ability  to  model  the  coupling  between  arc  dynamics  and 
processes  occurring  outside  the  capillary. 
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Figure  11.  The  Average  Molecular  Weight  of  50-bar  Plasmas  as  Computed  by 
CEA  and  the  SE  Model. 


5.  Conclusions 


The  chemical  compositions  of  plasmas  created  from  decomposed  polyethylene  as  computed 
by  (1)  the  NASA-Lewis  computer  program  for  the  calculation  of  complex  chemical  equilibrium 
(CEA)  and  (2)  an  SE  model  with  a  limited  reaction  product  set  have  been  compared.  Results  are 
presented  for  temperatures  from  5,000  to  20,000  K  and  pressures  from  10  to  200  bar,  conditions 
that  are  relevant  to  modeling  an  ablating-capillary  arc  discharge  during  the  ignition  phase  of  an 
interior  ballistic  cycle.  The  study  was  undertaken  to  establish  whether  CEA  could  be  used  to 
generate  source  species  concentrations  for  a  CFD  model  being  developed  to  characterize  plasma- 
propellant  interactions.  The  results  indicate  that  even  though  CEA  does  not  account  for 
Coulomb  interactions,  its  estimates  of  chemical  composition  should  be  adequate  for  this 
application. 
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Appendix  A: 

Fill  and  Cutoff  Procedures  Used  to  Establish 
the  Energy  Level  Structure  for  Carbon 
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The  calculation  of  the  electronic  partition  function, 


2,„(r)  =  E^,,y.*exp 


^  —  E’‘  ^ 
kj 

V  *  y 


(Al) 


for  plasma  specie  (i,f)  at  a  temperature  (T)  requires  knowledge  of  its  energy  levels  (Efj),  their 

statistical  weights  and  a  procedure  to  address  the  fact  that  the  sum  is  not  intrinsically 

bound.  In  a  plasma,  the  latter  concern  is  typically  resolved  by  tmncating  the  summation  at  a 
level,  max,  set  by 


<E~  -AE~ 

t.J  t.j 


(A2) 


where  E"j  is  the  ionization  potential  of  the  specie  and  A  E“j  is  the  reduction  in  the  ionization 

potential  resulting  from  its  interaction  with  charged  particles.  This  is  referred  to  as  a  cutoff 
procedure  and  will  depend  on  one’s  choice  of  ionization  potential  reduction  theory. 

Energy  level  information  for  atomic  species  can  be  found  in  a  number  of  references,  with 
Moore’s  tables^  being  perhaps  the  most  widely  used.  However,  usually  only  levels  with  low 
values  of  orbital  angular  momentum  are  known,  and  such  tables  are  almost  always  incomplete 
below  the  limit  given  by  equation  (A2).  In  the  case  of  (neutral)  carbon,  Moore’s  table  does  not 
include  any  information  on  f  or  higher  levels  of  orbital  angular  momentum,  and  information  on 
any  n  level  becomes  sparser  as  n  increases.  This  turns  out  to  be  a  significant  issue  for 
temperatures  above  about  15,000  K.  The  issue  arises  because  carbon  has  a  relatively  low 
ionization  potential — Ic,i  =  90,878  cm'^  (11.3  eV)  vs.  Ih,i  =  109,679  cm'^  (13.6  eV),  - 
117,345  cm'^  (14.5  eV),  and  Io,i  =  109,837  cm  '  (13.6  eV) — and  high  principle  quantum 
numbers  (with  statistical  weights  proportional  to  n^)  are  present  at  relatively  low  energies.'  (For 
example,  n  =  4  levels  of  carbon  are  observed  between  9.7  and  10.5  eV,  while  in  hydrogen,  these 

*Moore,  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol  1,  'H  to  ^^V,”  NBS- 
NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 
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levels  are  observed  at  12.7  eV.)  Since  these  “missing”  levels  have  high  statistical  weights,  their 
omission  can  lead  to  serious  errors. 

For  temperatures  above  about  15,000  K,  we  found  that  if  partition  functions  were  computed 
based  on  using  the  energy  levels  listed  by  Moore  alone,  they  were  significantly  lower  than  those 
tabulated  by  Drawin  and  Felenbok^  unless  AE“j  was  greater  than  about  4,000  cm'^  (0.5  eV).  As 

a  first  step  toward  completing  (or  “filling”)  Moore’s  table  for  carbon,  we  relied  on  the 
spectroscopic  data  tabulated  by  Striganov  and  Sventitskii  to  identify  “missing”  levels.  Then: 
assignments  of  the  data  in  their  tables  provided  the  4f  levels  (84,050  cm'^)  and  a  number  of  n  = 
5,  6,  7,  and  8  levels.  For  remaining  levels  (up  to  n  =  10),  we  estimated  energy  level  values  based 
on  trends  in  the  observed  levels.  For  example,  we  specified  the  energy  level  of  the  5g  levels  to 
be  86,570  cm"'  based  on  our  expectation  that  these  levels  would  be  higher  in  energy  than  the  5f 
levels  (86,470  cm'^),  but  by  a  value  less  than  the  separation  between  5f  and  5d  levels  (170  cm'^). 

Table  A-1  shows  the  data  used  for  computing  partition  functions  in  this  study.  (In  order  to 
reduce  the  size  of  the  table,  spin-multiplet  terms  from  Moore’s  table  are  combined  and  term 
values  limited  to  three  significant  digits,  but  five  significant  digits  were  used  in  the  actual 
computations.)  All  predicted  levels  with  energies  up  to  89,800  cm'^  are  included.  The  partition 
functions  generated  with  this  table  were  found  to  be  less  than  3%  different  from  those  tabulated 
by  Drawin  and  Felenbok^  for  all  and  temperatures  from  2,750  K  to  20,500  K.  The  highest 

energy  level  included  in  the  table  (89,800  cm"*)  corresponds  to  an  implicit  minimum  of 

1,100  cm"*  (0.15  eV).  If  n  =  11  levels  are  included — lid  being  observed  at  89,968  cm"* — 
they  begin  to  make  a  significant  contributions  to  the  partition  function  as  the  temperature 
approaches  20,000  K.  But  we  expect  that  such  levels  will  not  “exist”  for  plasmas  at  the  densities 
generated  by  an  ablating-capillary  arc. 

^Drawin,  and  P.  Felenbok.  Data  for  Plasmas  in  Local  Thermodynamic  Equilibrium.  Paris:  Gauthier- 

Villars,  1965. 

^Stiganov,  A.  R.,  and  N.  S.  Sventitskii.  Tables  of  Spectral  Lines  of  Neutral  and  Ionized  Atoms.  New  York: 
m/Plenum,  1968. 
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In  the  cases  of  and  the  levels  listed  in  Moore’s  table  proved  sufficient  for  the 
temperature  range  of  interest.  The  levels  employed  for  the  calculations  of  these  species  are 
shown  in  Tables  A-2  and  A-3.  The  decision  on  the  highest  level  to  include  for  this  study  was 
based,  in  part,  on  considering  that  only  levels  to  89,800  cm'^  were  retained  for  carbon. 
Comparison  with  Drawin  and  Felenbok’s^  data  provide  a  basis  for  the  adequacy  of  these  data 
sets  for  temperatures  to  20,000  K. 

In  the  case  of  hydrogen,  we  calculated  its  energy  levels  (in  wavenumbers)  and  their  statistical 
weight  using 


f 

=  109,679 

\ 


/ 


and 


8  HXk 


ir 


respectively.  Levels  up  to  ^  =  4  were  allowed.  Again,  comparison  with  the  partition  functions 
tabulated  by  Drawin  and  Felenbok  provided  a  basis  for  the  adequacy  of  this  approach  for 
temperatures  from  5,000  to  20,000  K. 


^Drawin,  H.-W.,  and  P.  Felenbok.  Data  for  Plasmas  in  Local  Thermodynamic  Equilibrium.  Paris:  Gauthier- 
Villars,  1965. 
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Table  A-1.  Energy  Levels  and  Statistical  Weights  Employed  to  Compute 
Partition  Functions  for  Carbon 


Design 

g 

1 

3 

5 

msm 

5 

msm\ 

1 

msm\ 

3 

mmsm. 

1,3,5 

3s  ‘P 

3 

0.0 


16.4 


43.5 


10,194 


21,648 


33,735 


60,300 


61,982 


64,100 


68858 


69,700 


70,744 


71,400 


72,611 


73,976 


75,300 


77,680 


78,100 


78,200 


78,300 


78,338 


78,531 


78,728 


79,300 


80,200 


80,564 


81,106 


81,300 


81,770 


82,252 


83,500 


83,730 


83,761 


83,800 


83,882 


83,949 


84,032 


84,050 


84,100 


84,852 


85,000 


85,000 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Stiganov  1968 


Moore  1948 


Moore  1948 


Moore  1948 


Author  Estimate 


m 


85,400 


85,626 


86,187 


86,300 


86,310 


86,400 


86,400 


86,400 


86,470 


86,500 


86,500 


86,650 


87,000 


87,632 


87,700 


87,800 


87,795 


87,800 


87,800 


87.807 


87,800 


87,831 


87,850 


8 


88,200 


88,500 


88,607 


Source 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Author  Estimate 


Author  Estimate 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 


Moore  1948 
Moore  1948 


Moore  1948 


Author  Estimate 


Author  Estimate 


Author  Estimate 


Moore  1948 


Moore  1948 


Moore  1948 


mnm 

3 

88,632 

Moore  1948 

mmm 

5,3,1 

88,639 

Moore  1948 

8s  * 

12 

88,560 

7d  'D 

5 

88,600 

Author  Estimate 

8p  * 

36 

88.850 

Author  Estimate 

7f-i  * 

480 

89,000 

Author  Estimate 

00 

720 

89,100 

Moore  1948, 
Author  Estimate 

9s  * 

12 

89,100 

Author  Estimate 

9P  * 

36 

89,300 

Moore  1948 

9d-k  * 

924 

89,500 

Moore  1948, 

4~VULUV/1  J-tJtlJJ  nnv.- 

10s  * 

12 

89,500 

Author  Estimate 

36 

89,650 

Author  Estimate 

lOd-l  * 

1152 

89,800 

Moore  1948, 
Author  Estimate 
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Table  A-2.  Term  Values  Used  for  Partition  Function  Computations* 


Designation 

g 

Level 

(cm'^) 

2p¥ 

2 

0.0 

4 

64.0 

2pP 

2 

43,000 

4 

43,022 

6 

43,051 

2p^D 

6 

74,931 

4 

74,933 

2pS 

2 

96,494 

"Moore,  C.  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol.  I,  'H  to 
^^V.”  NBS-NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 


Table  A-3.  Term  Values  Used  for  Partition  Function  Computations* 


Designation 

g 

Level 

(cm'^) 

1 

0.0 

2pP 

1 

52,315 

3 

52,338 

5 

52,395 

2p 

1 

102,351 

“Moore,  C.  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol.  I,  'H  to 
NBS-NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 
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Appendix  B: 

Thermodynamic  Properties  for 
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The  NASA-Lewis  computer  program  for  the  calculation  of  complex  chemical  equilibrium 
compositions'  (CEA)  establishes  the  equilibrium  state  of  a  system  by  minimizing  the  Gibb’s  free 
energy  of  a  reaction  product  set.  To  calculate  the  Gibb’s  free  energy  for  the  set,  it  relies  on  a 
library  of  thermodynamic  property  data  (THERMO.INP)  for  computing  the  heat  capacity 
[Cp(T)],  enthalpy  [H°(T)],  and  entropy  [^(T)]  of  each  product  at  a  temperature  (T).  The  library 
consists  of  a  set  of  product-specific  coefficients  (a,-,  that  yield  these  properties  via  the 
empirical  equations, 


C°(T) 

=  +a^T-' +a^+aJ  +  a,T^  +aj^  +aj\ 

R 


(Bl) 


H^iT) 

RT 


=  -aJ 


-2 


.^-1,  T  b, 

-|-  InT  +  o-,  +  a.  — ha, - h  a^ - h  a-, - 1 — , 

2  ^"2^3  4  '  5  T 


(B2) 


and 


S°iT)  T-^  r"  T* 

- =  — Cj - h  O2T  -h  03  InT  -h  a^T  -h  Oj - h  Cg - h  d-j - 1"  ^2  ’  (B3) 

R  2  2  3  4 

where  R  is  the  universal  gas  constant.  This  appendix  describes  the  determination  of  such 
coefficients  for  Though  concentrations  of  in  any  system  at  temperatures  <  20,000  K  are 
likely  to  be  very  small,  this  exercise  was  undertaken  in  the  interest  of  conducting  as  complete  a 
comparison  as  possible  between  CEA  and  the  SE  model  employed  by  Powell  and  Zielinski.^ 


'McBride,  B.  J.,  and  S.  Gordon.  “Computer  Program  for  Calculation  of  Complex  Chemical  Equilibrium 
Compositions  and  Applications  H.  User’s  Manual  and  Program  Description.’’  NASA  Reference  Publication  131 1, 
NASA  Lewis  Research  Center,  Cleveland,  OH,  1996. 

^Powell,  J.  D.,  and  A.  E.  Zielinski.  “Theory  and  Experiment  for  an  Ablating-Capillary  Discharge  and 
Application  to  Electrothermal-Chemical  Guns.”  BRL-TR-3355,  U.S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD,  1992. 
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As  a  first  step  toward  determining  the  (a„  coefficients  for  C  ,  Cp(T),  H^(T),  and  S^(T) 

3 

were  calculated  for  temperatures  from  298.15  to  20,000  K  using  the  standard  definitions, 


C;(T)  ^  ^.d^InZjT)  ^  ,^^lnZ(r) 

R  2  dT^  dT 


H‘’(T)  ^5  ^  ^dlnZ(T)  ^  H‘’(0K) 
RT  2  dT  RT 


(B5) 


(B6) 


S^jT) 

R 


1  +  lnl 

2 


kj 


27tk^T 

Ny 


s3/2\ 


+  T 


d]nZ(T) 

dT 


+  lnZ(r), 


(B7) 


where  Na  is  Avogadro’s  number  and  is  standard  pressure.  The  enthalpy  of  formation  for 
at  0  K  [/f’fO  AT)  =  4,150  kJ/mole]  was  determined  from  the  heat  of  formation  of  given  in  the 
NIST/JANAF  thermochemical  tables^  and  the  ionization  potential  of  given  by  Moore."^ 
Following  from  the  definition  of  the  partition  function,  the  terms  with  derivatives  were  computed 
per 


j.^lnZ(r) 

dT 


1 


kJZ{T)  k 


Igj£:*exp 


A 


KT 


(B8) 


and 


^Chase,  M.  W.  (ed.)  NIST-JANAF  Thermochemical  Tables.  4‘*'  edition.  Journal  of  Physical  and  Chemical 
Reference  Data,  Monograph  No.  9,  New  York;  American  Institute  of  Physics,  1998. 

“^Moore,  C.  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol  I,  *H  to  ^^V.” 
NBS-NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 
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,^dMnZ(T) 

dT" 


1 


k,TZ(T) 


kbT 


.-2  — 


1 


k,TZ(T)ir 


ESicE’'  exp 


k^T 


g^E  exp 


k.T 


yj 


.  (B9) 


The  energy  level  information  (£■*  and  g*)  needed  for  these  computations  were  extracted  from 
Moore’s  tables.'^  Table  B-1  shows  the  energy  levels  employed  for  the  computations.  Partition 
functions  computed  with  these  data  compare  favorably  with  those  tabulated  by  Drawin  and 
Felenbok.^  The  program  to  compute  thermodynamic  properties  was  checked  by  running  it  with 
the  data  listed  in  Appendix  A  for  carbon  and  verifying  that  the  thermodynamic  property  data  for 
carbon  tabulated  in  NIST/JANAF  thermochemical  tables^  were  reproduced. 


To  determine  the  a,  coefficients,  a  nonlinear  least  squares  fitting  routine  was  employed  to  fit 
equation  (Bl)  to  the  plot  of  Cp(T)IR  vs.  T  computed  via  equation  (B4).  The  a,  coefficients 
determined  in  this  manner  were  then  used  in  fits  of  (B2)  and  (B3)  to  (T)IRT  vs.  T  and  (T)/R 
vs.  T,  respectively — the  fits  yielding  bj  and  bz,  respectively.  Consistent  with  the  practice  in  the 
data  file,  coefficient  sets  were  separately  determined  for  the  temperature  ranges  298.15-1,000  K, 
1,000-6,000  K,  and  6,000-20,000  K.  The  coefficients  we  inserted  into  THERMO.INP  for  this 
study  are  given  in  Table  B-2.  We  note,  however,  that  this  table  is  primarily  for  informational 
purposes.  Gordon  and  McBride  use  routines  that  simultaneously  (rather  than  sequentially)  fit  the 
thermodynamic  property  functions.  Furthermore,  the  fits  are  constrained  so  that  the  polynomials 
coefficients  employed  on  either  side  of  a  transition  point  (1,000  and  6,000  K)  yield  the  same 
value  at  the  transition  point.  The  coefficients  are  also  determined  such  that,  within  certain  limits, 
the  polynomials  may  be  extrapolated  to  beyond  the  range  for  which  they  were  established.  We 
did  not  attempt  to  include  such  considerations. 


^Chase,  M.  W.  (ed.)  NIST-JANAF  Thermochemical  Tables.  4*  edition,  Journal  of  Physical  and  Chemical 
Reference  Data,  Monograph  No.  9,  New  York:  American  Institute  of  Physics,  1998. 

"^Moore,  C.  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol  I,  'H  to  ^^V.” 
NBS-NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 

^Drawin,  H.-W.,  and  P.  Felenbok.  Data  for  Plasmas  in  Local  Thermodynamic  Equilibrium.  Paris:  Gauthier- 
Villars,  1965. 
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Table  B-1.  Term  Values  and  Statistical  Weights  Used  in  Computation 
of  Thermodynamic  Property  Data* 


Designation 

g 

Level 

(cm-') 

28^ 

1 

0.0 

2pF 

1 

52315 

3 

52338 

5 

52395 

2p'P 

3 

102351 

“Moore,  C.  E.  “Atomic  Energy  Levels  as  Derived  from  the  Analyses  of  Optical  Spectra.  Vol.  I,  ’H  to 
^^V.”  NBS-NSRDS  35,  National  Bureau  of  Standards,  Washington,  DC,  1948. 


Table  B-2.  Coefficients  for  the  Computation  of  the  Thermodynamic  Properties  of  C^^ 


Coefficients 

298.15- 1,000  K 

1,000  -  6,000  K 

6,000  -  20,000  K 

ai 

2.5 

2.52243896E+04 

1.2021 1947E+08 

32 

0.0 

-7.51998881E+01 

-9.04691 343E+04 

as 

0.0 

2.58833405E+00 

2.91778474E+01 

34 

0.0 

-5.23213535E-05 

-3.85733714E-03 

35 

0.0 

1.65224996E-08 

2.75797361E-07 

36 

0.0 

-2.65014198E-12 

-8.555905 14E- 12 

37 

0.0 

1.697303213E-16 

9.89664769E-17 

bi 

4.991 70 12E+05 

4.9964775  lE+05 

1.19820991E+06 

b2 

2.5645969 

1.93672086 

-2.24190892E-K)2  | 
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